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Abstract 

We derive the nested Bethe Ansatz solution of the fully packed 0(n) loop 
model on the honeycomb lattice. From this solution we derive the bulk free 
energy per site along with the central charge and geometric scaling dimen- 
sions describing the critical behaviour. In the n = limit we obtain the 
exact compact exponents 7=1 and v = 1/2 for Hamiltonian walks, along 
with the exact value = SVS/A for the connective constant (entropy). Al- 
though having sets of scaling dimensions in common, our results indicate that 
Hamiltonian walks on the honeycomb and Manhattan lattices lie in different 

universality classes. 
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The configurational statistics of polymer chains have long been modelled by self-avoiding 
walks. In the low-temperature limit the enumeration of a single self-attracting polymer in 
dilute solution reduces to that of compact self-avoiding walks. A closely related problem isi 
that of Hamiltonian walks in which the self-avoiding walk visits each site of a given lattice 
and thus completely fills the available space. Hamiltonian walks are directly related to the 
Gibbs-DiMarzio theory for the glass transition of polymer melts [0. More than thirty years 
ago now Kasteleyn obtained the exact number of Hamiltonian walks on the Manhattan 
oriented square lattice 0. More recently this work has been significantly extended to yield 
exactly solved models of polymer melts 0] . The critical behaviour of Hamiltonian walks on 
the Manhattan lattice has also been obtained from the Q = limit of the Q-state Potts 
model [Q. In particular this Hamiltonian walk problem has been shown [0,^ to lie in the 
same universality class as dense self-avoiding walks, which follow from the n = limit in 
the low-temperature or densely packed phase of the honeycomb 0(n) model 

As exact results for Hamiltonian walks are confined to the Manhattan lattice, the be- 
haviour of Hamiltonian walks on non-oriented lattices, and the precise scaling of compact 
two-dimensional polymers, remains unclear [j^-Q. The exact value of the (Hamiltonian) 
geometric exponent 7^ was conjectured to be 7^ = ■y^, where 7^ = yI extracted via 
the Coulmb gas method for dense self-avoiding walks However, recent numerical in- 

vestigations of the collapsed and compact problems are more suggestive of the value 7^ = 1 

MM- 



More recently, Blote and Nienhuis [|T^ have argued that a universality class different to 
dense walks governs the 0(n) model in the zero temperature limit (the fully packed loop 
model). Based on numerical evidence obtained via finite-size scaling and transfer matrix 
techniques, along with a graphical mapping at = 1, they argued that the model lies in a 
new universality class characterized by the superposition of a low-temperature 0{n) phase 
and a solid-on-solid model at a temperature independent of n. This model is identical to 
the Hamiltonian walk problem in the limit n = 0. In this Letter we present exact results 
for Hamiltonian walks on the honeycomb lattice from an exact solution of this fully packed 



loop model. We derive the physical quantities which characterize Hamiltonian walks on the 
honeycomb lattice. These include a closed form expression for the connectivity, or entropy, 
and an exact infinite set of geometric scaling dimensions which include a conjectured value 



by Blote and Nienhuis ||T2[. Our results settle the abovementioned controversy in favour of 
the universal value 7^'^ = 1. 

In general the partition function of the 0(n) loop model can be written as 

^0(n)=E^^"^^^^^ (1) 

where the sum is over all configurations of closed and nonintersecting loops covering Mb 
bonds of the honeycomb lattice and M is the total number of lattice sites (vertices). Here 
the variable t plays the role of the 0(n) temperature, n is the fugacity of a closed loop and 
Ml is the total number of loops in a given configuration. 
For the particular choice t = tc, where 

tl = 2± 72^, (2) 

the related vertex model is exactly solvable with a Bethe Ansatz type solution for both 
periodic ||T3|-|r^ and open boundary conditions. This critical line is depicted as a 



function of n in Fig. 1. Here we extend the exact solution curve along the line t = 0, where 
the only nonzero contributions in the partition sum (1) are for configurations in which each 
lattice site is visited by a loop, i.e. with M = Mb ^7^. This is the fully packed model recently 



investigated by Blote and Nienhuis |12 



We consider a lattice of M = 2MN sites as depicted in Fig. 2, i.e. with periodic 
boundaries across a finite strip of width A^. The allowed arrow configurations and the 
corresponding weights of the related vertex model are shown in Fig. 3. Here the parameter 
n = s + = 2 cos A. In Fig. 2 we also show a seam to ensure that loops which wrap 
around the strip pick up the correct weight n in the partition function. The corresponding 
weights along the seam are also given in Fig. 3 [|TB|. We find that the eigenvalues of the 
row-to-row transfer matrix of the vertex model are given by 



J-i, sinh(e„ + if ) smh0^ sinh0^ 



where the roots 6^ and 0^ follow from 



e'^ 



sinh(6'Q, — 



n N 



'^^ sinh(^„ -<j)^ + i|) -p, sinh(^„ - - iA) 



[ sinh(^« + if)J ^i^i sinh(^„-(/.^-if) j-i,sinh(^„-^^ + iA)' 

(4) 

ie -pr sinh(0^ - - i|) ^ _ fr sinh(0^ - (pty - iA) ^ 
' iJ sinh(0, - ^« + i|) 3i sinh(0, - 0. + iA) ' ^ 

Here the seam parameter e = A for the largest sector and e = otherwise. Apart from 

the seam, this exact solution on the honeycomb lattice follows from earlier work by Baxter 

on the colourings of the hexagonal lattice |T9|. Baxter derived the Bethe Ansatz solution 

and evaluated the bulk partition function per site in the region n > 2. The corresponding 

vertex model was later considered in the region n < 2 with regard to the polymer melting 

transition at = ||2CI|| . 

More generally, the fully packed loop model can be seen to follow from the honeycomb 

limit of the solvable square lattice ^2^'* loop model Equivalently, the related vertex 



model on the honeycomb lattice is obtained in the appropriate limit of the ^2^^ vertex model 
on the square lattice in the ferromagnetic regime. This latter model is the 5^(3) vertex model 
||23| . One can verify that the above results follow from the honeycomb limit of the Algebraic 



Bethe Ansatz solution of the su{3) model with appropriate seam. It should be noted 



that Reshetikhin has performed similar calculations to those presented here, although 

in the absence of the seam, which plays a crucial role in the underlying critical behaviour. 

Defining the finite-size free energy as Jn = logAo, we derive the bulk value to be 

f°° sinh^Ax sinhfvr — X)x , , ^ 

/oo = / ^ — -^dx. (6) 

J-oo X sinhvrx sinhSAx 

This result is valid in the region < A < 7r/2, where the Bethe Ansatz roots defining the 

largest eigenvalue Aq are all real. We note that the most natural system size is a multiple 

of 3, for which the largest eigenvalue occurs with ri = 2A^/3 and r2 = A^/3 roots. In the 

limit A — s> /oo reduces to the known n = 2 value ||19|,[T3|, 



There is however, a cusp in the free energy at A = 7r/2. For A > 7r/2 the largest eigenvalue has 
roots 9a shifted by i7r/2. The result for /oo is that obtained from (^) under the interchange 
A ^ TT — A, reflecting a symmetry between the regions —2 < n < and < n < 2. Thus 
the value (0) holds also at n = —2, in agreement with the tc ^ limiting value |T3]. 



As our interest here lies primarily in the point n = (A = 7r/2), we confine our attention 
to the region < n < 2. At = 0, we find that the above result for /oo can be evaluated 
exactly to give the partition sum per site, as 

= 3^3/4, (8) 

and thus k, = 1.13975 . . . follows as the exact value for the entropy or connective constant 
of Hamiltonian walks on the honeycomb lattice. This numerical value has been obtained 
previously via the same route in terms of an infinite sum |2^. Our exact result (8) is to be 
compared with the open self-avoiding walk, for which fx^ = 2 + ^/2 and so /i = 1.84775 . . . 
It follows that for self-avoiding walks on the honeycomb lattice the entropy loss per step due 
to compactness, relative to the freedom of open configurations, is exactly given by 



- log 
2 ^ 



3^3 



-0.483161... (9) 



4(2 -Fx/2) _ 

The central charge c and scaling dimensions Xj defining the critical behaviour of the 
model follow from the dominant finite-size corrections to the transfer matrix eigenvalues 
| 26| . For the central charge. 



/^ - /oo + (10) 

The scaling dimensions are related to the inverse correlation lengths via 

= log(Ao/A,) ^ 2nCX,/N. (11) 
Here C = VS/2 is a lattice-dependent scale factor. 



The derivation of the dominant finite-size corrections via the Bethe Ansatz solution of 



the vertex model follows that given for the su{3) model in the antiferromagnetic regime E? 



(see, also [Q). The derivation is straightforward though tedious and we omit the details. 
In the absence of the seam, we find that the central charge is c = 2 with scaling dimensions 
X = AW + A(-), where 

A(±) = lgn''Cn+^ (h^fC-'h^ -]n h^, (12) 
8 o g 4 

C is the su{3) Cartan matrix and n = (ni,n2) with rii and n2 related to the number of 
Bethe Ansatz roots via ri = 2A^/3 — rii and r2 = — n2 The remaining parameters 

= {hf, hf) define the number of holes in the root distribution in the usual way [^. We 
have further defined the variable g = 1 — X/tt. 

With the introduction of the seam e = A, we find that the central charge of the fully 
packed 0(n) model is exactly given by 

c = 2-6{l-g)yg. (13) 



This is the identification made by Blote and Nienhuis [|T^. At n = we have g = 1/2, and 
thus c = — 1. On the other hand, both Hamiltonian walks on the Manhattan lattice [^,^ 
and dense self-avoiding walks lie in a different universality class with c = —2. However, 
as we shall see below, they do share common sets of scaling dimensions and thus critical 
exponents. This sharing of exponents between the fully packed and densely packed loop 
models has already been anticipated by Blote and Nienhuis in their identification of the 
leading thermal and magnetic exponents |T^. Here we derive an exact infinite set of scaling 



dimensions. 

Of most interest is the so-called watermelon correlator, which measures the geometric 
correlation between L nonintersecting self-avoiding walks tied together at their extremities 
X and y. It has a critical algebraic decay, 

{Mx)My))c-^\x-y\-^''\ (14) 
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where Xl is the scahng dimension of the conformal source operator 4>l{x) 0. As along the 
hne t = tc, these scahng dimensions are associated with the largest eigenvalue in each sector 
of the transfer matrix. The pertinent scaling dimensions follow from the more general result 

X=^-g{nl + nl-n,n,)-^-^^. (15) 

The sectors of the transfer matrix are labelled by the Bethe Ansatz roots via L = rii + n2. 
The minimum scaling dimension in a given sector are given hj rii = n2 = k for L = 2k and 
111 = k — l,n2 = k or ni = k,n2 = k — 1 for L = 2k — 1. Thus we have the set of geometric 
scaling dimensions corresponding to the operators (pL for the loop model, 

X2t-i = i J - t + l) - ^i^, (16) 

where k = 1,2, . . . The magnetic scaling dimension is given by X^ = Xi which agrees with 



the identification made in |T^. The eigenvalue related to X2 appears in the = 2 sector 
of the loop model, i.e. with two dangling bonds [|T^. At n = this more general set of 
dimensions reduces to 

X2k-i = \{k^-k), (18) 
X2k =\{k^~l). (19) 

In comparison, the scaling dimensions for dense self-avoiding walks are 
Thus we have the relations 

X2, = X™. (22) 

Note that Xi = X2 = and > for L > 2. Identifying X^ = X2 as in then the 
exponents 7 = 1 and u = 1/2 follow in the usual way These are indeed the exponents 
to be expected for compact or space filling two-dimensional polymers. 



The corresponding scaling dimensions for Hamiltonian walks on the Manhattan lattice 
are as given in (19) Q]. Exact Bethe Ansatz results on this model indicate that the scaling 
dimensions X^r = = 0, from which one can also deduce that 7 = 1 and u = 1/2 ||31|| . 
We also expect these results to hold for Hamiltonian walks on the square lattice. Extending 
the finite-size scaling analysis of the correlation lengths for self-avoiding walks on the square 
lattice 1^ down to the zero-temperature limit t = 0, we see a clear convergence of the central 
charge and leading scaling dimension to the values c = — 1 and Xi = for even system sizes, 
with X2 = exactly. These results are the analog of the present study on the honeycomb 
lattice where N = 3k is most natural in terms of the Bethe Ansatz solution. 

Our results indicate that fully packed self-avoiding walks on the honeycomb lattice have 
the same degree of "solvability" as self-avoiding walks on the honeycomb lattice. The fully 
packed loop model with open boundaries is also exactly solvable [^. The derivation of the 



surface critical behaviour of Hamiltonian walks is currently in progress. 
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FIGURES 

FIG. 1. The exact solution curve t = tc (2) of the 0(n) model on the honeycomb lattice as a 
function of n. In this Letter we extend the exact solution curve along the line t = 0. 

FIG. 2. The periodic honeycomb lattice of width L. Dashed lines indicate the position of the 
seam. 

FIG. 3. (a) The allowed arrow configurations and their corresponding vertex weights, (b) The 
modified vertex weights along the seam. 
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